FIELD-INDUCED SPIN EXCITATIONS IN 
RASHBA-DRESSELHAUS TWO-DIMENSIONAL ELECTRON 
SYSTEMS PROBED BY SURFACE ACOUSTIC WAVES 

P. Kleinert 

Paul-Drude-Intitut fur Festkorperelektronik, 
Hausvogteiplatz 5-7, 10117 Berlin, Germany 

V.V. Bryksin 
A.F. Ioffe Physical Technical Institute, 
Politekhnicheskaya 26, 194021 St. Petersburg, Russia 
(Dated: February 2, 2008) 

Abstract 

A spin-rotation symmetry in spin-orbit coupled two-dimensional electron systems gives rise to a 
long-lived spin excitation that is robust against short-range impurity scattering. The influence of 
a constant in-plane electric field on this persistent spin helix is studied. To probe the field-induced 
eigen-modes of the spin-charge coupled system, a surface acoustic wave is exploited that provides 
the wave-vector for resonant excitation. The approach takes advantage of methods worked out 
in the field of space-charge waves. Sharp resonances in the field dependence of the in-plane and 
out-of-plane magnetization are identified. 

PACS numbers: 72.25.Dc, 72.25.Rb, 72.10.-d 
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I. INTRODUCTION 



The perspectives of exploiting electron spin for information processing in spintronic de- 
vices stimulated basic research of spin-related phenomena.- In order to take advantage of 
traditional semiconductor technologies, the exclusive application of electric fields to gener- 
ate and manipulate a nonequilibrium spin density is of particular interest. In this respect, 
the spin-orbit interaction (SOI) provides an attractive mechanism for the pure electric ma- 
nipulation of moving spins in nonmagnetic semiconductors. Unfortunately, the very same 
SOI depends on the carrier momentum, which is randomized by elastic and inelastic scat- 
tering. This momentum dependence of the coupling leads to the dominant spin relaxation 
and dephasing in III-V semiconductors known as Dyakonov-Perel relaxation mechanism.- In 
general, however, details of spin dephasing depend on the involved scattering mechanisms, 
the band structure, and the crystal orientation (see Ref.- and references therein). 

An sizable anisotropy of the spin relaxation has been theoretically predicted^ for semi- 
conductor heterostructures grown along the [001] axis. According to these studies, a large 
spin-relaxation anisotropy should exist for the in-plane spin polarization. Moreover, un- 
der the condition that the strengths of Rashba (a) and Dresselhaus (j3) spin-orbit coupling 
constants become equal (a = (3), the anisotropy should be maximal. Based on this obser- 
vation a proposal for a novel spin-transistor device^ was put forward. Recent experimental 
data provided strong evidence for a huge in-plane anisotropy of the spin dephasing time in 
[001] GaAs/Al^Gai-^As quantum wells.-^ The efficient suppression of relaxation for spins 
aligned along the [110] orientation in heterostructures with equal Rashba and Dresselhaus 
SOI strengths is an important result that certainly stimulates further progress in spintronics. 

Other promising structures for spintronic applications are quantum well samples grown 
along the [110] direction. Very long spin- relaxation times on the order of nanosec- 
onds were experimentally demonstrated for the out-of plane spin polarization in such 
heterostructures.— 1 ^ This effect is caused by the absence of Dyakonov-Perel spin relaxation 
in the particular crystal direction. 

The common physical origin for weak spin relaxation in [110] quantum wells as well as 
in [001] heterostructures with a = /3 is the presence of an spin-rotation symmetry in the 
two-dimensional electron gas (2DEG). 12 This symmetry is robust and leads to an infinite 
spin lifetime and a persistent spin helix for an idealized model. The theoretical explanation, 
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which has recently been confirmed by experiment,— will initiate a number of further studies 
in this field. It is the aim of this paper to contribute to this development by treating 
weakly damped eigen-excitations of the coupled spin-charge system subject to a constant 
electric field. Our approach bears some similarities with a traditional field in solid state 
physics namely the study of space-charge waves in crystals. 14 The electric-field-mediated 
excitations of the coupled spin-charge system are probed by a simulated experimental set 
up, which provides the wave vector to resonantly excite internal eigen-modes of the biased 
SOI coupled 2DEG. We focus on excitation mechanisms that exploit surface acoustic waves. 

II. BASIC THEORY 

Our model is based on a Hamiltonian that describes conduction band electrons in a III-V 
semiconductor quantum well grown in the [001] direction, which is used to be the z axis of 
the coordinate system 

H = — h hw k ■ (T. (1) 

2m 

The effective electron mass, the in-plane wave vector, and the Pauli matrices are denoted by 
m, k, and Oi (i = x,y,z), respectively. The precession frequencies cjfc are due to the SOI, 
which contain both linear Rashba and Dresselhaus contributions 

wf= -(fc^-fc*), w£ = j(k x ,-k y ), (2) 

with a, (3 being the respective spin-orbit coupling constants. These vectors are also inter- 
preted as a momentum-dependent internal magnetic field that acts on the carrier motion. 
The neglected cubic terms to the frequency vector uj^ play only a secondary role in our 
study. Throughout this work, it is assumed that the spin-orbit coupling energy is much 
smaller than the Fermi energy Bp- The model is complemented by the inclusion of elastic 
scattering on a short-range impurity potential characterized by the elastic scattering time 
r. 

To elucidate the peculiarities of the spin-orbit coupled 2DEG, the coordinate system is 
rotated so that the [110] direction becomes a coordinate axis (k± = (k x ± k y )/y/2 — > r± = 
(r x ± r y ) I y/2). Along the spatial constant electric field E is applied. In addition, 

a surface acoustic wave (SAW) with the wave number Ksaw and frequency Q = v sav/Ksaw 
{vsaw is the sound velocity) should propagate in the same direction giving rise to an internal 
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electric field. The total electric field E + is expressed by 



E + (r + , t) = E + E S aw cos(K SAW r+ - fit) + 5E(r+, t), 



(3) 



is self-consistently calculated by taking into account the Poisson equation 

aifl(r + ,«) = 4g 
dr + e 

with F and n denoting the nonequilibrium and initially homogeneous carrier density, re- 
spectively (the components of the spin-density matrix F^, (fi = 1,2,3,4) are defined as in 
Ref.— ). e is the dielectric constant. Due to the excitation conditions in Eqs. ([3]) and (jlj), all 
quantities are independent of the spatial coordinate r_ for an infinite 2DEG. Along the r + 
axis, periodic boundary conditions are assumed. 

All electric field components act both on the charge and spin degrees of freedom leading 
to specific field-induced eigen-modes of the spin-charge-coupled system. The appearance 
and excitation of these modes is treated by drift-diffusion equations that are obtained from 
a kinetic theory for the spin-density matrix. For the given orientation of the constant electric 
field Eq and the SAW field, the set of equations decouples and takes the forrn^ 



dt drl 



dE + F h TJr dF_ 

F + ^-2^+^7 
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2 

2' 
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dE+F- 2 F. 
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T s - nE 



E + - 
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2D d 



dE F 

F + + a + + - K + I uE + - 2D 
dr. 



H dr 
d 



dr. 



(5) 

(6) 
(7) 

dr + y~ dr ' * ' 

where F± and F z denote the in-plane and out-of plane nonequilibrium spin densities, respec- 
tively. D is the diffusion coefficient, fi the carrier mobility, and K + = 2m(a + (3)/h 2 a wave 
number constructed from the spin-orbit coupling constants. In addition, we introduced the 
spin-relaxation times 



Fz + ^ d _E±El + K+ (^ E+ -2D- d 



CI. 



- = 4DK 2 , — 



cos 2 <p 



1 



sin 2 ip 



(9) 



with K = m\J a 2 + (3 2 /h . The angle cp is used to quantify the coupling constants according 
to a/^/a 2 + j3 2 = cos(y? + 7r/4) and j3/ \f a 2 + j3 2 = sm(ip + vr/4). For the limiting cases 



4 



of a pure Rashba and Dresselhaus model, we have ip = — 7r/4 and 7r/4, respectively. When 

the SAW field is switched off, the applied constant electric field E induces a homogeneous 

in-plane magnetization F_ that is calculated from 

/q) /-cos 2 <p ^ T ^ „ dn , , 

sm ip dEF 

Taking into account the Einstein relation \x = eDdhm/dep, this magnetization is expressed 
by the result = HKfiEodn/ 'dsF (with K = ma/h 2 ) for a Rashba model (/3 = 0) that 
was derived by Edelstein.— When the SOI coupling strengths become equal (a = j3), the 
spin-relaxation time t s _ diverges so that the result for F^ strongly depends on the initial 
condition and/or additional spin relaxation mechanisms quantified by the scattering time 
r sp (t~} — > r~} + T^ 1 ). The divergency is also circumvented in the ac response.— 1 ^ In 
all cases, the field-induced in-plane magnetization disappears for a = (3. In Eq. (JSJ), G z 
describes the uniform generation of an out-of plane spin polarization by optical means or by 
the application of a perpendicular constant magnetic field. 

From Eqs. (jl]) and (jSJ), an expression for the total time- dependent charge current density 
I(t) is obtained 

e d5E ,, . . ^ ^ ^ dF ,-fiK ^ , . 

j(r+,t) + = I{t), j(r+,t) = efxE + F - eD— - V2 cos<^eF_, 11 

47r at or + m 

which completes the set of basic equations necessary to study electric-field-mediated effects 
in the spin-charge coupled 2DEG. The drift-diffusion Eqs. (jSJ) to (jHJ) decouple into two sets 
of equations for the components F, F_ and F z , F + . Accordingly, the spin-charge coupling 
prescribed by Eqs. ([SD and (jB]) is treated in the next Section, whereas Eqs. ([7D and (JSJ) are 
solved in Section 4. 



III. SPIN-CHARGE COUPLING 

The solution of Eqs. (jSJ) to (jSJ) is facilitated by the fact that all quantities depend only 

onz = Ksaw t + ~ Qt. The derivative with respect to z is denoted by a prime. Introducing 

the dimensionless functions / = F/n, f± = F±/n, and f z = F z /n and the parameters 

v E + v 6E o v *r 1 4vre DK SAW 

y+ = ~et, Y = — , Vl E = fiE K SA w, — = fin, A = — , (12) 

we obtain the following set of coupled ordinary differential equations 

- Qf_ - AQ E fl + — /_ + n E (Y + f_)' + — f {0) (Y+f - 2Af ) = 0, (13) 



-fir M r + r + /-Af-v / 2-^rCos^_ = -, Q e t m Y' = f - 1. (14) 

These equations are solved by applying perturbation theory with respect to the SAW electric 
field Esaw- When the SAW field is absent (Esaw = 0), all quantities are independent of z 
and we obtain for the constant charge current density 



jo 1 



y/2 - — cos ipf^ 



(15) 



m/iEo 

with jo = enfiE Q . The spin-induced contribution on the right-hand side of this equation 
results from the homogeneous in-plane magnetization that vanishes for a system with 
equal Rashba and Dresselhaus coupling constants (a = (3). 

To proceed, we calculate the stationary current contribution 51 associated with the SAW 
electric field. This time-independent current contribution arises from carriers that are driven 
by the SAW. A similar current can also be generated by a moving optical lattice in a 
photorefractive crystal. As a periodic boundary condition along the r + axis is assumed, 
the solution is searched for in the form of a Fourier series f(z) = n + Sf(z), Sf(z) = 
^ p exp(ipz)<5/(p). Taking into account the Fourier transformed version of Eq. fllll) 



Jo 

we obtain 



6 P , = - ip nr M Y(p) + y+(P ~ P')f(p) - ipA/(p) - V2 



hK 

mfiEo 



cosp/_(p), (16) 



61 

jo 



Y f 



SAW 



l-v^-^cos^' 



[6f(l) + Sf(-1)) 



(17) 



mnE 

with Ysaw — Esaw/Eo. The nonequilibrium fluctuation of the charge density 5/(1) 
for p — 1 is expressed by the variation of the internal electric field SY(1) via 6Y(1) = 
—i6f(l)/(Q E TM)- The latter quantity is calculated from the Fourier transformed Eqs. f|T3|) 
and dE]), which are given by 
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—ip{VL — Vt E + ipAQ E ) + 



6f-(p) 



(18) 



(0) 
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-— f m Y S Aw(S p ,i + S p ,- 
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hK 
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(19) 



Even at the absence of SOI, the SAW field alone gives rise to a charge current that is easily 
calculated from Eqs. (fT8l), (fT9]h and (jT7j) 

( 5I\ Yg AW Vt E Vt — Vt E 



_ = isaw_Ei _E: (20) 

This current contribution is due to space-charge waves, the dispersion relation of which is 
obtained from the poles in Eq. (|20|) 

n = fiE K SAW + dk Iaw^) > (21) 

where tm denotes the Maxwellian relaxation time. The physical origin of this mode, which 
has been thoroughly studied in the field of photorefractive materials^, are oscillations of 
the free electron gas. In the same way as the current density jo has a spin complement 
in Eq. ( fl5i) . the well-known SAW-induced charge current in Eq. ( 120]) has a related spin 
contribution, which is also calculated from Eqs. (TT71) and ({IS]) . (jl9l) with the result 
fSl\ h K + K SAW f { ^Yl AW 



\jo J 2 4m 



2 1 

SAW) 




2-2 1 



+ — (22) 

T s - T M J 

2 \ /n n n2 



The SAW-induced electronic current in Eq. fT2U|) is independent of the spin-degree of freedom 
and disappears at the resonance frequency Q = Q E of space-charge waves. This result has to 
be compared with Eq. (122]) that holds for the spin complement of the SAW-induced charge 
current, which may exhibit a sharp resonance at Q = Q E . Besides space-charge excitations 
that also appear in Eq. (120]) . there exists a new eigen-mode, in which r^ 1 is replaced by 
the spin-scattering rate 2/r s _ that vanishes for a = (3. This spin-mediated mode has the 
character of a space-charge wave with an infinite Maxwellian relaxation time (for a = (3). 
In experiment, the SAW related stationary charge currents in Eqs. ( 1201) and ( 122]) can be 
distinguished by their qualitative different dependence on the electric field E . 



IV. SPIN-SPIN COUPLING 



For the considered set-up not all four components of the spin-density matrix couple to each 
other. According to the special field orientation, the coupling between the spin components 
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F z and F + is separated from the charge density F and the in-plane spin component F_ . The 
related drift-diffusion Eqs. (J7J) and (jSJ) take the dimensionless form 

- nf + - + A /+ + n E ( Y+ UY - ^Y + f z + 2*£**±f K = 0, (23) 

T s+ saw ft SAW 

- nf z - m E f« + A/, + + - 2^M/; = ^, (24) 

r s J<saw Ksaw n 

where G z denotes the source of an out-of-plane spin polarization. Without the application 

of the SAW field, we obtain the steady-state solution 

f (o) _ Gz l n (Q) _ nE Q (0) 

Jz 2/T s + (fiE y/D , J+ DK + Jz ' 1 J 

The out-of-plane spin polarization /i "* is suppressed by the electric field.— By determining 
this suppression in experiment, the spin-relaxation time t s can be determined in a way 
similar to the well established experimental set up that is based on the Hanle effect (c.f., 
for instance, Ref.-). Furthermore, the electric field Eq induces an in-plane spin polarization 

/f. 

As in Section 3, the effect of the SAW is treated by a perturbation expansion f z = fz°^+5f z 
and /+ = + Sf + . The Fourier transformed basic Eqs. (l23|) and ( |24l) are easily solved, 
when non-linear perturbations are neglected. The solution for the lowest-order Fourier 
components is given in the Appendix. 

To determine the steady-state spin polarization under the mutual influence of a constant 
electric field and the SAW field, we focus on the most interesting case a = (3 and assume 
Qe t m 1. Under these conditions, the solution given by Eqs. (|3T|) to (133]) simplifies and 
takes the form 

*/+(0) . Y* AW Sl E i n E (l-q 2 )Z 1 + n(l-Al-2iqA + ) \ 

" 2(i + A^) ne \ (n-n + )(n-n_) /' 1 ' 

Sfz(0) Yhw^E j n E (l-q 2 )Z 2 + i Q(l-Al + 2 i qA) \ 

f jp) 2A(i + A 2 + ) Ke | y 1 } 

with the shorthand notations Z\ = 3A+ - 1 + zA(l - A+), Z 2 = A (3 - 2iA + A 2 + ), q = 
K + /K S awi an d A + = qA. The eigen-frequencies of the spin excitations are expressed by 

fi± = vE (K SAW ± K + ) - iD(K SAW ± K + f. (28) 
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100 200 300 400 

E [V/cm] 

FIG. 1: SAW- induced stationary in-plane spin 
polarization 5f x (0) referred to fx^YsAW for 
^SAW = 5-6 /im, [i = 2.3 10 4 cm 2 /Vs, and D = 
24 cm 2 /s. Results are shown for K + /Ksaw = 
0.8, 0.85, and 0.9. 



100 200 300 400 

E a [V/cm] 

FIG. 2: SAW-induced stationary out-of-plane 
spin polarization Sf z (0) referred to ff^YsAW 
for Xsaw = 5.6 /im, fi = 2.310 cm 2 /Vs, 
and D = 24 cm 2 /s. Results are shown for 
K + /K S aw = 0.9 and 0.925. 



The appearance of resonances at O = Q± that occur in the in-plane and out-of-plane field- 
induced spin polarization is the most remarkable result of this paper. Due to the spin- 
rotation symmetry, a soft mode develops in the system that can be probed by a SAW, 
which provides the wave vector Ksaw for the resonant excitation. The low-frequency mode 
Q — > becomes increasingly undamped for Ksaw —> ±K±. As a remnant of these coherent 
long-lived spin excitations, a weakly damped resonance appears, when the applied constant 
electric field E satisfies the condition fl = v saw Ksaw — ^E (K S aw ~ K + ). 

Numerical examples are shown in Figs. [T] and [2] for the SAW-induced in-plane and out- 
of-plane spin polarization. The parameters used in the calculation refer to GaAs quantum 
wells.— ^ It is to be noted that in the displayed relative spin polarizations, the electric field 
Eq enters via the field- dependent quantities and /J ^ calculated from Eq. (I25j) . Both 
components of the magnetization exhibit a sharp resonance in the electric field dependence at 
fiE = vsaw/ (1 — K+l Ksaw) as indicated by the vertical line in Fig. [21 With increasing ratio 
K + / Ksawi the resonance is shifted to higher field strengths and becomes more pronounced. 
Under the condition Esaw ~ Eq and at the resonance, the out-of-plane spin polarization 
may significantly exceed the spin generation given in Eq. (1251) . As the suppression of 
at high electric fields can be compensated by an appropriate in-plane magnetic field 2 ^, 
the experimental demonstration of the effect should be feasible. Its most salient feature is 
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the strong variation of the field-induced magnetization in the vicinity of the resonance field 
strength. This switching of the magnetization caused by an applied constant electric field 
could be useful for future spintronic device applications. 

V. SUMMARY 

Recently long-range transport of spins has been experimentally demonstrated in [001] 
GaAs quantum wells with balanced Rashba and Dresselhaus terms as well as in Dresselhaus 
[110] quantum wells. The idealized model that describes both systems exhibits an exact spin 
rotation symmetry with an associated soft mode and a persistent spin helix. This long-lived 
spin coherence state is transferred to a field-dependent eigenmode of the system, when a 
constant in-plane electric field is applied. To resonantly probe the field-induced excitation, 
the required wave- vector must be provided by an appropriate experimental set up. Similar to 
the well established physics of space-charge waves one can produce optical gratings for that 
purpose. In this paper, we treated an excitation mechanism via an acoustic wave propagating 
along a [001] GaAs quantum well, in which both Rashba and Dresselhaus SOI exist. Due 
to the spin-charge coupling, specific eigenmodes are identified both in the charge current 
density and the in-plane and out-of-plane spin polarizations. The SAW-induced change of 
the stationary charge current is composed of two contributions. The first one is independent 
of the SOI and exhibits a pole that is due to oscillations of the free carrier density. In the 
second contribution, which is related to the field-induced homogeneous spin accumulation, 
a new eigenmode appears with the dispersion relation Q = fiE Q K S Aw ~ (DKg AW + 2/r s _). 
In this equation, the spin-scattering time r s _ (which diverges in the case a = f3) takes over 
the role of the Maxwellian relaxation time tm that is responsible for the damping of pure 
charge oscillations. 

Characteristic field-dependent eigenmodes appear also in the SAW-induced components of 
the magnetization. Most interesting is the case, when the Rashba and Dresselhaus coupling 
strengths become equal. According to the dispersion relation Q = a undamped soft mode 
occurs, when the wave-vector K S aw of the SAW approaches the shifting vector K + = 2m(a+ 
(3) I h 2 . As a remnant of this persistent spin helix, a sharp resonance occurs in the electric-field 
dependence of the in-plane and out-of-plane magnetization at /iE = vsaw/(1 — K+/Ksaw)- 
In the vicinity of this resonance, the spin polarization drastically changes with a slight 
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variation of the electric field strength. It is a challenge to experimentally verify this field- 
induced switching of magnetization in a 2DEG with balanced Rashba and Dresselhaus SOI. 

Solution of kinetic equations In this Section, the coupled Eqs. ([23]) and ffZ$§ are solved by 
a perturbation approach with respect to the SAW electric field Ysaw- First, the equations 
are expressed in terms of a Fourier series, which leads to the result 



-i P n + P 2 An E + — + i P n E ] sf+(p) + tpYsA ^ E [5f + (p + 1) + sf+(p - 1)] 

T s+ J 2 

A '+" £ 1 - 2i P A)S M p) - *+»W"' lSW + 1) + SUP - U] 

ZK 



K 



SAW 



SAW 



SAW 

2% 



(S p ,i - 5 p -t) - ipY(p) 



+ 



K + f° z n E 



K 



SAW 



SAW 



(S P ,i + (J Pj _i) + Y{p) 



(29) 



-ipn + p 2 aq e + - + i P n E } 8f z (p) + ipYsA ^ nE [ 5 f z {p + i) + 8f z ( P - 1)] 

Ts J 2 

- 2ipA)5f + (p) + K f* AwQE [5 f + {p + 1) + 5f + (p - 1)] 
j^saw ^j^saw 

Ysaw 



(30) 



mi 



2% 



-(<W - Sp-i) - ipY(p) 



Y f 



SAW 



{S P ,i + <5 Pj _i) - Y[p) 



Ksaw 

To calculate the steady-state solution, the equations for the p = components 5f+(0) and 
<5/z(0) are solved. We obtain 

Ysaw 



Sf+(0) 
Sf z (0) = 



l + 2(« 2 + /3 2 )AV(« + /?) 2 
Ysaw 



Re[A + 5f z (l)-8f + (l)\, 
Re [A+*/+(l) + */,(!)]. 



(31) 



(32) 



l + 2(« 2 + /? 2 )AV(a + /5) 2 
In a second step, the Fourier components Sf + (1) and 5/ z (l) have to be calculated. Disre- 
garding nonlinear contributions, the remaining equations are easily solved with the result 

ci A s - c 2 B cj. c 2 A s+ + ciB 



A, 



n = I Y(l) 
2 



XsAW 



aa+ + s 2 



(33) 



Ysaw 



)( 



Ksaw J 



n 



+ - iA - — 



A s+ =-?— + z(l-zA--^-), £ 

"E7"s+ -ft SAW 



(1 - 2zA). 



Vt E r s Q E 

The SAW-induced stationary magnetization of the biased 2DEG is expressed by Eqs. (|3l 
to 
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